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ABSTRACT 

We present Atacama Large Millimeter/submillimeter Array (ALMA) 860 /im imaging of four high- 
redshift (z=2.8-5.7) dusty sources that were detected using the South Pole Telescope (SPT) at 1.4 mm 
and are not seen in existing radio to far-infrared catalogs. At 1.5^' resolution, the ALMA data reveal 
multiple images of each submillimeter source, separated by 1-3'', consistent with strong lensing by 
intervening galaxies visible in near-IR imaging of these sources. We describe a gravitational lens 
modeling procedure that operates on the measured visibilities and incorporates self-calibration-like 
antenna phase corrections as part of the model optimization, which we use to interpret the source 
structure. Lens models indicate that SPT0346-52, located at z=5.7, is one of the most luminous and 
intensely star-forming sources in the universe with a lensing corrected FIR luminosity of 3.7 x IO^^Lq 
and star formation surface density of 4200 M© yr~^ kpc~^. We find magnification factors of 5 to 
22, with lens Einstein radii of 1.1—2.0'' and Einstein enclosed masses of 1.6-7.2x10^^ M©. These 
observations confirm the lensing origin of these objects, allow us to measure the their intrinsic sizes 
and luminosities, and demonstrate the important role that ALMA will play in the interpretation of 
lensed submillimeter sources. 

Subject headings: Gravitational lensing: strong — Galaxies: high-redshift — Techniques: interfero- 
metric 
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1. INTRODUCTION 



Half of the energy produced by all objects in the his- 
tory of t he universe has been absorbed and reemitted 
by dust ( [Dole et al. 2006). The Cosmic Infrared Back- 
ground, lirst de tected by the Cosmic Background Ex- 
plorer satellite ([Pu^et et "aL 1996 Hauser et al. 1998 



Fixsen et al.|| 199*8 1, is the aggregate emission tr om indi- 



vidual dust y galaxies across cosmic time (e.g., Lagache 
et al. 2005|). The brightest of these dusty star-forming 



galaxies (DSFGs) were discove red in deep submihimeter 
wavelength images of the sky ([Smail et al.||1997 Hughes 
et al. 1998[ [Barger et al. 199op , and have lunimosities 
in excess of 10^^ L(?) emitted primarily at rest wave- 



length in the far-infrared. With star formation rates 
> 100 - 1000 M0 yr-\ this population of DSFGs con- 
tributes a significant fraction of the total star formation 
density of the uni verse at z ^ ^ ~ 3, wh ere their abun- 
dance peaks (e.g., 'Chapman e t al.||200' 5^. These objects 
are the progenitors of the massive galaxies we observe 
today. 

Despite the enormous total luminosity of the brightest 
DSFGs, their detection at submihimeter wavelengths re- 
quires lengthy exposures for ground-based facilities, and 
they are generally quite dim at optical/NIR wavelengths 
due to extinction. Studies of these objects and their 
extreme star formation rates are limited by the obser- 
vational costs of observing all but the brightest spectral 
lines and the poor spatial resolution achievable compared 
to the typical size of the star-forming regions. Gravita- 
tional lensing provides a solution to both of these prob- 
lems, as has been demonstrated in a few spectacular cases 
(e.g., Kn eib et al . 2004; Swinbank et al. 2010; Riechers 
pt al |20lf| " l^u et al.^20 12). Lensed starburst galaxies can 
be examined at higli spatial resolution and with a more 
diverse s et of diagnostics t han the unmagnified popula- 
tion (Swi nbank et al.||2010| ). 

Predictions of a large population of gravitationally 
lensed, hig h redshift DSFGs (Blain 1996; Negrello 



et al.||2007) were recently v erifiedTpy large-area millime- 



ter/ subnnllimeter surveys (Vieira et al. 20101 Negrello 



eTa^ [Wardlow etldp01:j| ). [Hezaveh Holder 

(|2011 j also predicted the number counts of bright lensed 



objects for mm- wavelength surveys using a detailed nu- 
merical method, with a proper treatment of finite source 
effect and lens ellipticities, confirming that realistic lens 
models were able to match the observed number counts 



of dusty sources reported in Vieira et al. (2010). These 



galaxies have a sky density of ^ 0.1 deg ^, and therefore 
can only be found in large numbers in extensive surveys. 
The South Pole Telescope (SPT; Carlstrom et al.,2011) , 
which surveyed 2500 square degrees to ^ mJy depth at 
wavelengths of 3, 2, and 1.4 mm, has provided a sample 



of about one hundred candidate lensed sources (Vieira 



et 

foum 



ID 

ina 



2010). Initial investigations of these objects have 
them to have properties consistent with unlensed 
starbursts, except for their large apparent luminosities 
(Greve et al. 2012). Morphological evidence of lensing 
cannot be discerned in dat a from the SPT su rvey or the 
single-aperture followup of Greve et al. (2012), except in 
rare cases of lensing by clusters ol galaxies, so arcsecond- 
resolution submihimeter imaging is required. 

Paris-Sud XI & CNRS, 91405 Orsay Cedex, France 



The A tacama Large M illimeter /submihimeter Array 
(ALMA, |Hills et al.|2010[ ) has begun operation in Chile, 
providing unprecedented submihimeter sensitivity even 
in early science. In this work, we employ ALMA to mea- 
sure the arcsecond-scale structure of dusty extragalactic 
SPT sources at millimeter wavelengths, confirming the 
lensed nature of the four sources presented here. The 
observations reported here represent <10% of our Cy- 
cle sample and use only the compact configuration 
data, which was delivered first. Nevertheless, from these 
ALMA data we are able to model the lensing geometry 
of these sources and de-magnify them, allowing them to 
be placed in the proper context within the high-redshift 
galaxy population. We are also able to infer the total 
mass and ehipticity of the lenses, a first step toward 
using the lensed submihimeter emission to character- 
ize the lensing potential and its substructure. In Sec- 
tion 2 we describe the ALMA observations and support- 
ing data and in Section 3 we describe a modeling tech- 
nique for interferometric measurements of gravitationally 
lensed sources. Additional details on the integrated self- 
calibration step are included in the Appendix. In Section 
4 we discuss the properties of the sources and lenses and 
present the conclusions in Section 5. Throughout this 
work we assume a ACDM cosmology, with WMAP7 pa- 
rameters, with h = .71, Qm = 0.27, and Qa = 0.73 
dKomatsu eFaI]|2011 ). 



2. OBSERVATIONS 
2.1. ALMA Imaging 

The primary observations for this work were obtained 
from ALMA under a Cycle program (2011.0.00958.S; 
PI: D. Marrone) in which 47 sources identified in the 
SPT survey are each observed in both the compact and 
extended array configurations. The first data release for 
this program includes 20 sources observed in the compact 
array configuration, here we focus on the four sources for 
which these low-resolution data were sufficient to resolve 
the targets into multiple components. Of the remaining 
16 sources, at least 8 are not point-like at the resolution of 
these observations, but we defer lens modeling for these 
sources until the remaining data are in hand. 

The sources were targeted for brief snapshot observa- 
tions with the dual-polarization Band 7 (275-373 GHz) 
receivers on two dates, 2011 Nov 16 and 28. The first lo- 
cal oscillator was set to 343.8 GHz, with all four spectral 
windows configured in time domain mode with 128 chan- 
nels of 15.625 MHz width centered at 5.125 and 7 GHz IF 
in each sideband. There were 16 and 14 antennas avail- 
able on these days, respectively, arranged in a compact 
configuration. The total elapsed time (for all sources, 
including those not published here) in the observations 
was approximately 4.1h. The total integration time per 
source was 61 and 91 seconds in the first and second 
tracks, respectively, with 6.1-second sampling of the vis- 
ibility data. The array alternately observed the science 
targets and gain calibrators (30 seconds), observing the 
calibrator every 3—4 minutes. Additional sources with 
known positions, precisely established against the Inter- 
national Celestial R eference Fr ame using very long base- 
line interferometry (Ma et al. n.998), were added to the 
tracks to verify astrometry and calibration and observed 
with the same cycle as the science targets. 
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The flux scale was set with observations of Callisto on 
the first day. On the second day the flux scale was de- 
rived by setting the flux of quasar J0403-360 to 1.84 Jy, 
as reported by the ALMA staff from adjacent calibration 
observations. The absolute flux density scale is correct 
to within 15%. The antenna gains are equalized through 
gain calibration (amplitude and phase) on the main cal- 
ibrators in each track. Short-timescale phase correc- 
tion is achieved using the ALMA water vapor radiom- 
etry (WVR) system. Very little variation is observed in 
these gain amplitudes through the tracks, and there is 
no evidence for atmospheric decorrelation on the longest 
baselines in the phase scatter of the target or calibrator 
visibilities after WVR phase correction. The data were 
processed with th e Common Astronomy Software Appli- 
cations package (Mc Mullin et a l.|2QQ7[|Pe try et al.|2Q12D 
using standard steps for aTc^ntinuum observcition. 

The four well-resolved sources are listed in Table [l} 
we refer to them throughout the paper with shortened 
versions of their coordinate-based names (e.g., SPT0346- 
52 for SPT-S J034640-5205.1). Observations on Nov 16 
included baselines of 15-150 kA, resulting in synthesized 
beams of 1.5''xl.3'' (FWHM) for the RA=5^ sources. 
The uv coverage on Nov 28 was less uniform and spanned 
15-240 kA, providing a synthesized beam of 2.f xO.9'' for 
the other two sources. Deconvolved source images and 
beam shapes are shown in Figure [T] 

2.2. Redshift Determinations 

Table [T] shows source and lens redshift s, as available, 
for these four SPT sources. Redshifts for three of the 
dusty sources were obtained in another ALMA Cycle 
project (2011.0.00957.S; PL A. Weiss) through a spectral 
line survey conducted using the Band 3 (84-116 GHz) re- 
cei vers. Comp lete results of this survey will be reported 
in |Vieira et al.| ( |2012| ) and |Weifi et al.| ( |2012[ ). For each 
of these sources, multiple high-signiticance lines are de- 
tected in the Band 3 spectral scan, providing unambigu- 
ous redshifts. In the case of SPT0538-50, a redshift was 
determined from a combination of millimete r- wavelength 
and optical spectroscopy, as described in [Greve et al. 
( [2012) , and the source was not included in the ALMA 
redshift search proposal. 

The combination of NIR pre-imaging and submillime- 
ter interferometric observations described below pro- 
vided the basis for ground-based spectroscopic obser- 
vations of the putative lens galaxies. The fields of 
SPT0346-52, SPT0418-47, and SPT0529-54 were tar- 
geted first with i?-band pre-imaging, then with multi- 
object masks, using the Mask Exchange Unit of F0RS2 



(Appenzeller et al. 1998| ) on the Very Large Telescope 
(VLT) UTl, and exposmg them for 3x900 sec integra- 
tions each. The data were collected during February and 
March, 2012, in Service Mode (ESO program ID 088.A- 
0902) under an average seeing of ^ l.l'^ The galaxies 
were observed with slits of V' in width using the G300 V 
grism, yielding a velocity resolution of ~ 650 km s~^ or 
13A, sampled at - 3.3A pixer^ SPT0538-50 was ob 



served with the XSHOOTER echelle spectrograph (Ver- 
net et al.|2011 ) in longslit mode. For these observations, 
the slit widths were V' (UV-B), 0.9'' (VIS-R), and 0.9'' 
(NIR); the corresponding resolving powers are R = 5100, 
8800, and 5600, sampled with 3.2, 3.0, and 4.0 wave- 
length bins respectively (after on-chip binning by a fac- 



tor of two of the UV-B and VIS-R detectors). 6x900 sec 
integrations were obtained in February, 2010, however 
half of these were taken under worse conditions and we 
use only the better three integrations in the spectrum 
presented here. 

The observations were prepared and the data reduced 
using the standard ESO pipelinj^ performing bias and 
flat corrections, background subtraction, fringe correc- 
tion, registration and combination, wavelength calibra- 
tion, and Id spectral extractions. Spectra are shown in 
Figure (2] With the exception of SPT0346-52, redshifts 
for the Tensing galaxies were measured by fitting Gaus- 
sians to the Ca K+H absorption lines. 

2.3. Other Supporting Data 

The small primary beam size of the 12m ALMA an- 
tennas at 350 GHz required that we improve upon the 
~10" positional uncertainty of the SPT detections before 
proposing ALMA observations. These targets were ini- 
tially followed up with the Submillimeter Array (SMA; 
Blunden|20Q4l|HQ et aH2004 ), or the Australia Telescope 
Compact Arra}j^(ATCA; SPT0529-54 only). The SMA 
observations (Botnwell et al. in prep.) had typical angu- 
lar resolution of 3"xl2" because of the low declination 
of the sources, but those data gave some indication that 
the submillimeter emission was resolved. The 100 GHz 
ATCA detection of SPT0529-54 was of low significance 
and did not reveal any ringlike structure. 

Deep near-infrared (NIR) imaging data were acquired 
from the Southern Astrophysical Research (SOAR) Tele- 
scope (OSIRIS; Depoy et al. 1993), the Very L ar^e Tele- 



scope (VLT ISAAC; [Moorwood et al.||1998[ ), and the 



Hubble Space Telescope^ as part of lollowup programs to 
examine lens properties and identify rest-frame optical 
emission from the background sources. 

The luminosities in Table [T] are derived from photo- 
metric observations at millimeter to submillimeter wave- 
lengths. In addition to SPT photometry (1.4 and 2 mm 



we use 870 /im observations from LABOCA (Siringo 
et al. 1 12009) on the APEX telescope (in Max Planck 



time) and Herschel-SFIRE photometry at 250, 350, and 
500 /im. The LABOCA data were analyz ed according 



to the procedure described in Greve et al. (i2012}. The 

rwgiFi — 



Hersc hel data are reduced as described in Iweip et al.l 



(2012). All photometric data, including S PT measure- 



ments and the data previously reported in Greve et al. 
("2012^ for SPT0529-54 and SPT0538-50, are provided 
in WeiB et al. (2012). The LABOCA measurements of 
the total 870 /im flux density agree with the total flux 
density in the modeled ALMA data (described below) 
to within calibration uncertainties, typically 10%. Lumi- 
nosities are calculated using a g reybody model like that 
described in [Greve et al.| (|2012 ) , though in the present 
case the r = 1 wavelength is^ parameter of the fit. The 
effect of freeing this parameter is to broaden the peak 
of the spectral energy distribution (SED) as needed to 
match the short-wavelength Herschel data points. 

3. RESULTS AND ANALYSIS 

'^^ http://www.eso.org/sci/facilities/paranal/instruments/fors - 
VLT-MAN-ESO-19500-1771 

http:/ /www. narrabri.atnf.csiro.au/ 
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Fig. 1. — ALMA 350 GHz CLEANed images of the lensed DSFGs (contours), on top of NIR imaging of the galaxy field (greyscale) 
from Hubble Space Telescope, VLT, and SOAR. Left to right they are SPT0346-52 (HST) and SPT0418-47 (VLT) in the top row and 
SPT0529-54 (SOAR) and SPT0538-50 (HST) in the bottom row. The ALMA contours begin at 5cr, increasing in steps of lOcr in the top 
left and bottom right images, and 5a in the other two, as in Figure pi The synthesized beam is shown in the lower left of each image. The 
HST images are a composite of the FllOW and F160W filters, the VXT images are K band from ISAAC, and the SOAR image is K band 
from OSIRIS. 



Figure [TJpresents the resolved structure of the 350 GHz 
emission associated with the SPT sources, as revealed in 
Cycle ALMA observations, and the NIR emission at 
the same positions. In all cases there is a clear NIR 
counterpart to the submillimeter source, though with no 
structural correspondence between the infrared sources 
and the emission found in the ALMA images. The mor- 
phologies of these sources at submillimeter wavelengths 
are indicative of gravitational lensing. The redshift iden- 
tifications reported in Table [l] with very high redshifts 
measured for the dusty emission behind low-redshift fore- 
ground galaxies, clearly confirm these systems as gal axy- 
galaxy lenses. A similar finding was reported in Ne- 
grello et al. (2010*), where bright sources selected at much 
shorter wavelength from wide-field submillimeter surveys 
with Herschel-SFIRE were also confirmed to be gravita- 
tionally lensed. 



3.1. Lens Modeling 



The ALMA interferometer measures visibilities, 
Fourier components of the sky intensity distribution 
across a two-dimensional range of spatial frequencies, 
rather than directly imaging the emission. To prop- 
erly compare these data with a source model, we must 
perform our analysis in the visibility plane, where the 
measurement and its noise are well understood. Past 
tech niques for modeling in terferometric lens observations 
(e.g., Wucknitz et al. 2004) have generally operated on re- 
constructed images that are subject to difficult-to-model 
biases and noise properties. Furthermore, residual er- 
rors in calibration in the visibility data, such as those 
arising from imperfect knowledge of the antenna po- 
sitions or uncompensated atmospheric delay, are often 
corrected as part of the imaging proces s through an it- 
erativ e clean/self-calibration technique (Cornwell et al.| 
|1999|). However, the inclusion of the cleaning step m 
the determination of these corrections, which are then 
applied to the visibility data themselves, changes the 
data in ways that are not easily included in the mod- 
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TABLE 1 
Source and Lens Parameters 



















Source Intrinsic Properties 




ID 








Ml 






^1/2 


^FIR 


^FIR 


'S'l .4mm 


'S'seO/^m 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


SPT-S J034640-5205.1 


5.656 




1.124±0.004 


3.73±0.04^ 


0.55±0.01 


5.4±0.2 


0.59±0.03 


37.3±2.8 


24 


8.1 


23.0 


SPT-S J041840-4752.0 


4.224 


0.265 


1.390±0.012 


2.39±0.04 


0.20±0.03 


21.0±3.5 


1.07±0.17 


3.8±0.7 


0.74 


1.6 


5.0 


SPT-S J052903-5436.6 


3.369 


0.140 


1.536±0.017 


1.64±0.04 


0.10±0.03 


9.4±1.0 


2.39±0.24 


3.8±0.5 


0.15 


3.8 


15.6 


SPT-S J053816-5030.8 


2.782 


0.404 


1.987±0.009 


7.15±0.05 


0.13±0.02 


20.5±4.0b 




4.5±0.9 




1.5 


5.8 


SPT0538-50 A 












19.8±4.6 


0.52±0.12 


3.0^ 


2.4^ 






SPT0538-50 B 












21.9±3.7 


1.61±0.33 


1.4^ 


0.12^ 







Note. — Column 1: SPT source name. Column 2: Background source redshift. Column 3: Lens redshift. Column 4: Einstein radius (arcsec). 
Column 5: Lens mass, interior to te (lO-'^-'^M©). Column 6: Lens ellipticity. Column 7: Total magnification of background source. Column 8: 
Source radius, determined as the half-width at half maximum for the Gaussian source component in the model fit (kpc). Column 9: Intrinsic far 
infrared luminosity (IQ-^'^Lq). Column 10: Intrinsic source fiux (luminosity per area; 10-*^^ L0/kpc^). Column 11: Intrinsic 1.4 mm fiux density, 
obtained as the SPT fiux density divided by /i (mJy). Column 12: The same as (11), but for the ALMA 350 GHz (860 /im) fiux density. Parameter 
uncertainties do not include a contribution from the cosmological parameters. An additional 4% uncertainty in mass is found by marginalizing over 
the WMAP7 parameter Markov chains. 

^ Assuming that the lens is located at z = 0.8, see section^ 

^ Total magnification of the two source components, see section^ 

^ Derived assuming that L^ir is divided between components in the same ratio as the fiux density in the model. 
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Fig. 2. — Optical spectra of the foreground lenses, from the VLT. The positions of major lines are marked by thin dashed lines. Sky 
lines are marked with grey shaded regions. 



eling uncertainties. Here we describe a visibility-based 
lens modeling technique that simultaneously determines 
these self-calibration phases so that we incorporate the 
full range of uncertainty present in the measurements. 

We model the lenses by generating model lensed im- 
ages that are subjected to simulated observations and 
compared to the data. The source is assumed to have a 
symmetric Gaussian light profile with four free parame- 
ters: flux density F^, radius R^, and positional offsets 
from the lens center X^, Ys. The lens is modelled as a 



Singular Isothermal Ellipsoid (SIE) with five free param- 
eters: mass inside the Einstein radius Ml, ellipticity e^, 
orientation angle 6l (east of north), and position X^, Yl. 
The SIE profile has been shown to be a good approxima- 
tion to galaxy-scale density profile s (Treu & Koopmans 



20041 iKoopmans et al.||2006H2009'). ^ ince lensed imase 



positions provide very precise measurements of the pro- 
jected mass interior to the images we report this robustly 
measured quantity in Table |T] the total mass (Ml) inside 
the Einstein radius (r^;) of tne lensing galaxies. However, 
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the total halo masses associated with these galaxies will 
be much higher. Magnification, /i, is calculated as the 
ratio of the total lensed to unlensed flux. 

Given a set of lens and source parameters we make a 
high-resolution image of the lensed source, which we pad 
with zeroes and Fourier transform to generate model vis- 
ibilities. For each ALMA visibility, we interpolate the 
Fourier transformed image to the ALMA v coordi- 
nates. We also correct for the antenna primary beam 
attenuation by multiplying the sky model images by 
the primary beam pattern before sampling the Fourier 
modes. We use a symmetric Gaussian with FWHM of 
18^' for the primary beam, which leads to minimal atten- 
uation for these well-centered sources. 

The agreement between data and model visibilities is 
determined by calculating the between them, which 
requires an estimate of the visibility noise. Because we 
observe strong sources, the visibility scatter has a con- 
tribution from the sky signal and cannot be used to de- 
termine the noise level. Instead, we derive noise lev- 
els by measuring visibility scatter after differencing vis- 
ibilities that are adjacent in time for the same base- 
line/polarization/IF, which has the effect of removing 
all sky signal. This gives results that are identical to 
those found by scaling the visibility noise to obtain a re- 
duced of unity for the best-fit models. We explore 
the model parameter space using a Markov Chain Monte 
Carlo (MCMC) method with Metropolis-Hastings sam- 
pling. An example of the parameter degeneracies in these 
fits is shown in Figure |3] 

An additional complication of the modeling procedure 
is the presence of uncorrected antenna-based phase er- 
rors in the data. Observations of test sources, quasars 
with positions referenced to the International Celestial 
Reference Frame (ICRF), show significant residual phase 
errors after the primary phase calibration step. For sim- 
ple source structures, self-calibration using CLEAN com- 
ponents as input for the phase correction is a standard 
procedure to improve image fidelity (e.g., Taylor et al." 
|l999). However, in this application it would add signifi- 
cant complication to model the clean/self-cal process as 
part of the fitting, and leaving the phase errors uncor- 
rected can significantly bias the model parameters be- 
cause of the strong sensitivity to image flux ratios in the 
model. We have therefore developed a procedure to de- 
termine the self-cal phases as part of the model fitting. 
We optimize the fo^" each step in the Markov chain 
by adjusting the N — 1 antenna phases. We find that 
the resulting phases vary little over the chain and closely 
resemble those found for nearby point sources added to 
the tracks to test the calibration and astrometry, giving 
us confidence in this method. Additional details of the 
method and simulations of its effectiveness are provided 
in the Appendix. 

The model parameter space is complex, with the pos- 
sibility of multiple isolated minima separated by high 
barriers in x^. To decrease the possibility of missing 
important minima in the posterior we search the space 
more broadly by "tempering", which is similar to the 



simulated annealing method (e.g., [Press et al.||2QQ7 ). A 
control parameter, T (the analog ot temperature), is in- 
troduced to flatten the posterior surface and to make the 
minima more accessible. This is achieved by raising the 
posterior to the power 1/T (with T>1). 



4. DISCUSSION 

A primary goal of this work is to determine the lens- 
ing configuration for bright SPT starburst galaxies and 
derive the total magnification. Using a simple multi- 
comp onent source model a nd a range of lensing geome- 
tries, iHezavetTeFaL (2012) showed that differential mag- 
nification of a DSFG can distort the SED in unpre- 
dictable ways by differentially lensing different source 
plane regions (see also Serjeant 2012^ for a similar ef- 
fect concerning molecular line ratios). A lens model is 
therefore essential if we wish to correctly interpret ob- 
servations of these targets and place them in context 
with existing samples of unlensed starbursts. In prin- 
ciple, spatially resolved imaging of each molecular line is 
required to map that particular line to the source plane, 
however a single lens model based on continuum imag- 
ing, combined with physically motivated models for the 
relative filling factors of the emitting regions in other 
bands or molecular lines, can be used to place limits on 
the differential magnification between components. 

Models for the four sources are shown in the rightmost 
panels of Figure |4l Key parameters of the models are 
reported in Table^ The lens models permit the calcu- 
lation of some intrinsic properties of the lensed galaxies 
by simply dividing observed properties by the magnifica- 
tion. The four sources presented here were detected with 
1.4 mm flux densities ranging from 30—45 mJy. By cor- 
recting for the lensing magnification (Figure pi) we find 
that the intrinsic 1.4 mm and 860 /im fluxes of these 
sources vary by a factor of several. The 860 /im fluxes 
span the range for sour ces identified in b l ank field sur- 
veys with SCUBA (e.g.JGreve e t al."2004! 'Coppin et al^ 
|2006, Scott et al.|2008| [WeiJj et a l. 2009), suggesting tliat 
even without the lensing boost the galaxies presented in 
this work would still be identified as luminous starbursts. 
The far-IR luminosities are even more widely varying, 
though always ultra luminous infrared galaxy (ULIRG) 
class, indicating star formation rates of several hundred 
to several thousand Mq yr~^. 

The intrinsic sizes of the source-plane emission regions 
{Ri/2 in Table [ij the HWHM for the Gaussian model 
components) range from 0.5 to 2.4 kpc. This is at 
the lower range of the sizes generally inferred for star- 
bursts from a variety of observables (e.g., Bothwell et al. 



20101 J Rujopak arn et al.|[2QTT| |Daddi et al...2010, 1'l'acconi 



et al. 



2006, 201Q[ a ndTeferences therein), includings: syn- 



chrotron emi ssion ([Chapman et al. | |2004 Biggs & Ivison 
20081) low-J (llvison et al.||201l|) and even hi^her-J (e.R., 



Laccbni et al. 2006) CO measurements. However, the se- 
lection of strongly lensed galaxies can result in a biased 



distribution of intrinsic source size, as noted by [Hezaveh 



et al.| ( |2012| ). 
"Some comments on the individual sources: 
SPT0346-52: This source has the lowest magnification 
of the set (5.4±0.2), and being the brightest of the four 
targets in apparent 5'i.4mm, it is by far the most lumi- 
nous after correction for magnification. Unlensed, with 
a flux density of 26 mJy it would be among the 860 /im- 
brightest high-redshift galaxies known. Given the source 
size, the flux (luminosity per area) is2.4xl0i^ L0 kpc-^ 
implying a star- formation rate of 4200 M© yr~^ kpc~^ 
assu ming a standard conversion for starbursting galax- 
ies (Kennicutt 1998). This is remarkably high, 50 x 
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Fig. 3. — Parameter degeneracy plot for SPT0346-52, a representative example of the model uncertainties. Some nuisance parameters 
are not shown. Marginalized distributions for each parameter along the bottom axis are shown as histograms on the diagonal. The bottom 
row shows correlations between /i, which is derived from the model and not a fit parameter, and the model parameters. The contours show 
the 1, 2, and 3 a confidence regions. 

generally more compact than high-redshift galaxies of 
similar luminosity, suggestive of a self-regulating process. 
Because we find the luminosity of SPT0346-52 to arise in 
a very compact region, its mode of star formation may 
be more similar to lower-redshift ULIRGs than to most 
other high-redshift starbursts, despite being at z = 5.7. 
Walter et al. (2009) found a similar star formation den- 
sity over a region of similar size in a z=6A quasar, so this 
is not unprecedented in the early universe. The optical 
spectrum of the lens does not show any lines or features 
that result in a robust redshift for the lens. However, 



higher than the aver ag^e value found in m arginally re- 
solved starbursts by iTacconi et al. (20061) and a fac- 
tor of several higher than observed for individual GMCs 
in the highl y magnified source SMM J2135-0102 (Swin- 



bank et al. 2011). The flux is comparable to the Ed- 



dington limfTl I'hompson et al. 2005), and the degener- 
acy between /a and K^- (shown m l^lgure [3| is such that 
larger magnification decreases and further increases 
the flux. 



Thompson et al. (2005) noted that such fluxes are rel- 
atively common among low-redshift ULIRGs, which are 
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"17.50'* 17.00-^ 16.50^ 
RIGHT ASCENSION (J2000) 



Fig. 4. — Modeling of the observed submillimeter emission. In each row, the panels from left to right are: the "dirty" image, model dirty 
image, residuals after subtraction, CLEANed image, and fully resolved lens model. Contours in the left two columns start at ±3 times 
the rms noise in the residual map (third column) and increase in steps of 10 times the rms noise, except in the second row where they 
increase in steps of 5 times the rms noise. Contours in the residuals are zb 1, 2, etc., times the rms noise, which is (top to bottom) 1.5, 
1.7, 1.0, and 1.0 mJy. The contours in the CLEANed images start at 5 cr and increase in levels of lOcr for the top and bottom rows and 5cr 
in the two middle rows. The apparently high significance structure in the dirty maps away from the source is due to the sidelobes of the 
synthesized (or "dirty") beam, and should therefore be reproduced by the model in the second panel. The insets in the last panels show 
a magnified view of the positions of the emission in the source plane (greyscale) relative to the lensing caustics (red). In the source model 
for SPT0538-50 the greyscale is truncated at 10% of the peak intensity to make the second source component visible. The intensity scale 
for the right panel, in Jy arcsec"^, is given by the colorbar in each figure. 

the only lens modeling parameter that is degenerate with 
the redshifts is the lens mass. The mass of this lens is 
reported assuming z = 0.8 for the lens. 

SPTO4I8-47 Siud SPT0529-54: These sources should 
present nearly complete Einstein rings at higher res- 
olution. In the case of SPT0418-47, this is due to 
the excellent alignment of source and tangential caustic. 
SPT0529-54 is the most extended source in this sample 
and its large size compared to the caustic fills the ring 
effectively. In both cases, the luminosities, submillime- 
ter fiuxes, and source sizes are very comparable to those 
observed in other samples of distant starburst galaxies 



(e.g.,"Rujopakarn et al."2011). The posterior of the model 
parameters for SPT0529-54 show two separate peaks in 
the lens and source positions, but the magnifications of 
both models are similar. The extended array data on 
this source can possibly break this degeneracy. 

SPT0538-50: The large angular size of this ring sug- 
gests a massive lens, with the models indicating that the 
lens has a projected mass of nearly 10^^ within a 
radius of 10 kpc. Close alignment between the compact 
source (labeled component "A" in Table fl]), and caustic 
again leads to a large magnification. A single-component 
source does not provide an adequate fit, leaving a signifi- 
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to correct for the magnification of the sources presented 
in this work and derive intrinsic properties, finding the 
galaxies to be typical high-redshift, DSFGs. The sen- 
sitivity of ALMA permits these lens models to be con- 
strained in short observations. Longer observations of 
lensed starbursts in future cycles will therefore enable 
studies of ISM structure and lower luminosity molecular 
lines that are otherwise impossible to observe in unlensed 
systems. 
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Fig. 5. — The source magnification distributions derived from 
the lens models. In the case of SPT0538-50, the total magnification 
of the two source-plane components is shown. 

cant residual structure to the south east of the lens center 
(Figure [6|. A far better match to the data is the model 
shown in Figure [II (bottom), which includes a second 
source component (labeled component "B" in Table [T]), 
offset from and much more extended than the first, and 
representing 30% of the total source-plane luminosity. 

5. CONCLUSION 

We have used ALMA to image the submillimeter emis- 
sion from four DSFGs discovered by the SPT. We find 
that these objects resolve into ring-like structures ex- 
pected from gravitational lensing, a picture confirmed 
by the redshift information we have for the submillime- 
ter emission and NIR counterparts. We present a visibil- 
ity modeling procedure to fit gravitational lens models 
to these data and simultaneously correct the unknown 
phase errors of the antennas introduced by, e.g., imper- 
fect antenna positions. From this technique, we are able 

APPENDIX 
LENS MODEL SELF-CALIBRATION 

Interferometric phase calibration procedures generally leave some residual phase errors due to imperfect baseline 
solutions, uncompensated atmospheric delays, or other effects. The magnitude of these errors depends on many 
factors, such as calibration interval and calibrator-source separation, and their importance depends on the signal to 
noise ratio in the data and the complexity of the imaging task. In the present application, small phase errors can 
redistribute flux between lensed images and without properly accounting for such effects in our lens modeling, the 
derived model parameter distributions may be significantly in error. 

A standard procedure to correct antenna-based phase errors is self-calibration. Using a source model derived from 
images of the corrupted visibilities, phase corrections are derived, the source is re-imaged, and another iteration can 
be made using a new source model produced from the phase-corrected visibilities. For the present purposes, the most 
significant disadvantage of this phase correction scheme is that the uncertainties associated with the source model 
against which the data are self-calibrated are not included in the lens modeling, and any structural errors introduced 
by noise or phase errors may become a permanent part of the lens model. Here we propose a method to incorporate the 
phase correction into the lens modeling procedure, using the lensed structure as the source model for self calibration. 

To implement this self calibration technique, we use a perturbative approach in which we assume that the current 
data is equal to the model plus an antenna based phase corruption. We can write the expression for foi" the visibility 
phases between the data and model as 



X 



(Al) 



where is the phase of the z'th visibility, is the phase difference between the model and the data and S(j)k is the 
phase delay in the /c'th antenna. d^i/d(j)k is a Jacobian matrix containing the gradients of the observed visibilities 
with respect to changes in antenna phases. For N antennas and M visibilities d^/d(j) is an x M matrix whose z/c'th 
element is 1 if the first antenna of z'th visibility is A:, —1 if the second antenna of z'th visibility is k and zero otherwise. 
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5^38-17.50-^ 17.00'^ 16.50'^ 
RIGHT ASCENSION (J2000) 

Fig. 6. — The best fit lens model for SPT0538-50 using a single source component, for comparison with Figure Left to right are the 
"dirty" image, the model dirty image, the difference between data and model, and the full-resolution source modeL Contours in the two 
panels on the left start at ±2 times the rms noise in the residual map (cr), increasing in steps of lOcr, while contours in the third panel 
increase in steps of zblcr. The peak in the residual map has a significance of >5cr. As shown in Figure ^ a second source component 
completely eliminates this high peak and has a significance >10 a in the visibility space. 

C is a diagonal covariance matrix whose nonzero elements are approximately {o'xy/\V\)'^^ where axy is the rms error 
on the real or imaginary part of each visibility and \ V\ is the visibility a mplitude. Note that this approximation to the 
covariance matrix is only valid in the limit of high signal to noise data dWrobel fc Walk^pOO] ). 

To minimize we set its derivative with respect to antenna phases to zero {dx^ /d(j)i = 0). This allows us to write 
the antenna phase offsets as 



(A2) 



where F is the Fisher matrix calculated as 



(A3) 



At every iteration of the MCMC code, the value of x? is minimized for the postulated model by deriving calibration 
phases using equation 



A2 



The resulting x^ is used to evaluate the likelihood and determine the next link in the chain, 
thereby incorporating the uncertain phase correction in the parameter exploration. In the limit of an intrinsically 
Gaussian distribution for these phase parameters, this is equivalent to marginalizing over the phases. In Figure [7] we 
compare the results of this procedure with the standard self-calibration based on CLEAN components, and the raw 
data. The improvement in source subtraction is significant, even compared to the standard CLEAN procedure. 

A second test of the simultaneous fitting of the lensed emission and antenna phases is shown in Figure [8 Simulated 
observations of a typical lens model were created with realistic noise levels, and antenna phase errors added to these 
visibilities. The MCMC fitting algorithm was applied to these multiple realizations of the same source to verify that 
the antenna phase errors are recovered. The excellent agreement between input and recovered phase demonstrates 
that the simultaneous fitting of the lens model and antenna phases does not bias the antenna phase measurement, 
despite the complicated source structure of the lens models. 
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4^18^^40.5^ 40'^ 39.5"' 39. 

RIGHT ASCENSION (J2000) 

Fig. 7. — A comparison of image residuals for three different corrections to the antenna phases. The dirty image of SPT0418-47 is shown 
in greyscale beneath each panel. Contours show the image residuals after the source model is subtracted from the dirty image. Left: 
Residual structure after the model that fits the raw data (no self-calibration) best is subtracted from the data. Middle: The same, but 
for data calibrated using a CLEAN-based self-calibration. Right: The same, but for data calibrated using the procedure described in this 
work. The presence of residual structure away from the source, with no corresponding residual at the source position, is a clear sign of 
imperfectly corrected phase in the left two panels. 
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Fig. 8. — Antenna phase errors recovered in 70 simulated observations of a lens model with different noise realizations. Phase errors 
were introduced to the antennas, with values indicated by the vertical dashed lines. The histograms show the phase values recovered in fits 
to the simulated visibilities following the technique outlined here. 
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HIGH RESOLUTION OBSERVATIONS OF SPT0346-52 

The models in Figure [4] were derived using the compact configuration data that were available at the time of 
submission. Higher resolution data were delivered later, and permit a direct comparison of the model based on the 
low-resolution data with the higher resolution observations of SPT0346-52. Figure [9] shows the predictions of the best 
fit model (from fitting to the compact data, presented earlier in this work) for the uv-coverage of the extended data (red 
contours). The black contours show the extended data. The contours demonstrate a high degree of agreement between 
the predictions and the new observations. A lens model for the extended configuration data shows a consistent model, 
resulting in magnification (from modeling the extended data alone) of 5.26 ±0.12 in agreement with the magnification 
derived from the compact data. 




3M6"'41.30* 41.0* 

RIGHT ASCENSION (J2000) 

Fig. 9. — A comparison of the model proposed for SPT0346-52 from the low-resolution observations presented above and the higher- 
resolution observations obtained after submission. The "dirty" image of the extended configuration observations are shown in greyscale 
and black contours. The red contours are the predicted appearance given the model of Figure [4] and the uv sampling of the new data. 
Contours are drawn at -5 (dashed) and 20, 40, 60 (solid) times the rms noise. 
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